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We analyze the effects of sequences of standing wave pulses on a Bose-Einstein condensate (BEC). 
Experimental observations are in good agreement with a numerical simulation based on the band 
structure theory in the optical lattice. We also demonstrate that a coherent control method based 
on such sequences of pulses is very efficient for experimentally designing specific momentum states. 
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I. INTRODUCTION 

Atomic interferometry is very useful in fundamental 
studies of coherence, decoherence and phase shifts and 
for practical precision measurements, with the example 
of gravimeters, gyroscopes, and gradiometers [l|. Bose- 
Einstein Condensate (BEC) based atomic interferometry 
provides high contrasts, long integration times and the 
possible use of small devices [H, [3[. In such atomic in- 
terferometry, coherent momentum manipulation is very 
effective for splitting and recombining the condensate [J- 
Q, hence realizing the interference. In some precision 
measurements, the accumulated phase is positively corre- 
lated to the atomic velocity, so that the larger the atomic 
momentum is, the more precise the measurements can be 
if the measuring time is unchanged. 

The diffraction of atoms from standing wave light, 
which is usually divided into three regimes as Bragg, 
Raman-Nath and channeling regimes [7] by the interac- 
tion intensity and duration, is one of the common meth- 
ods to coherently manipulate the momentum states of 
a condensate 

[iHlii- One-pulse Bragg or Raman-Nath 
scattering can be applied for momentum states' prepara- 
tion, but those techniques are then limited by constraints 
on the pulse properties [7]. A momentum manipula- 
tion method by multi-pulse standing waves is proposed 
in where the momentum states can be designed, 

but not observed yet, and the pulses are still restricted 
in the Raman-Nath regime. 

In this paper, we apply a method for flexible manipula- 
tion of the atomic momentum states, where the standing 
wave pulses are less limited in pulse intensities and du- 
rations. The atomic diffraction from one, two, three and 
four standing wave pulses are demonstrated in our exper- 
iments and systematically analyzed by the band struc- 
ture theory of one-dimension optical lattice. With this 
method, we are able to design and realize several specific 
momentum states, which may be applied in atomic in- 



terferometry. In principle, this method could be used for 
designing a wide range of possible target states. 

This paper is organized as follows. In Section II, a the- 
ory to interpret the scattering process by a sequence of 
pulses is presented, where the standing wave is treated 
as a one-dimensional optical lattice. We derive a concise 
expression for calculating the probability of each momen- 
tum state at the end of the process. In Section III, the 
experiments of one, two, three and four pulses' scatter- 
ing are demonstrated and compared to theoretical sim- 
ulations. It is found that a correction due to momen- 
tum dispersion may be introduced into the theory for a 
better agreement with the experimental results. Section 
IV presents the experimental realization of several useful 
momentum states by coherent control. Section V con- 
tains the discussions and the conclusions. 



II. THEORETICAL MODEL 

We first consider a non-interacting condensate being 
diffracted by a sequence of square shaped standing wave 
pulses with the successive durations Ti{i = 1, 2, 5 + 1), 
separated by the intervals Tfi{i = 1, 2, s). The stand- 
ing wave consists of a pair of laser beams far-detuned 
enough to suppress the spontaneous emission. 

The periodic potential (one-dimension optical lat- 
tice) [isl, Q introduced by the ac Stark shift can be 
described as V{x) = [/q cos^(/clx), with the trap depth 
Uq and the laser's wave vector kL = 27r/AL (A^ is the 
wavelength of the laser). The lattice leads to a band 
structure of the energy spectra, of which the eigenvalues 
of the energy En^q and eigenvectors |n,g) (Bloch states) 
are labeled by the quasi-momentum q and the band index 
n, and they satisfy the equation: 



H\n,q) = En^q\n,q) 



(1) 
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where the Hamiltonian H = /2M + UQCOs^ikLx)^ with 
the atomic momentum p and the atomic mass M . The 
Bloch states form a quasi-momentum space. In the lat- 
tice, the spatial periodicity of the wave function results 
in separated peaks in momentum space, corresponding 
to the reciprocal lattice vector 2k l- 
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When a condensate with an initial momentum = 
h{q + 2mokL) {h is the Plank constant, — 1 < <7 < 1 , 
mo = — 1, 0, 1, ...) is abruptly loaded into a lattice, the 
wave packet can be described as a superposition of the 
Bloch states: 



l«'(^ = 0)) = £|n,(z)(n,9Ko) 



(2) 



where (n, q\pr^ 



(mo). The n Bloch state evolves 



independently as e '^^ri,qt/h^ total wave function 

evolves as 



iEn,qt/h\ 



n,q) 



(3) 



n=0 



While the incident light is switched off after the duration 
Ti, the wave function is projected back to the momentum 
space from the quasi- momentum space. The coefficient 
6(mo,m,ri) of each \prn) state (m = ...,—1,0,1,...) can 
be acquainted as: 



6(mo,m,ri) 



n=0 



Cn,(?(^o)Cn,(?(m)e' 



-iEn,qTl/h 



(4) 



For a zero initial momentum of the condensate, the sub- 
script q can be omitted for simplification and mo = 0. 
For one pulse scattering, the population of the \pm) state 
is = |6(0,m,ri)p. It can be seen that the probabil- 
ities of the momentum states after one scattering pulse 
depend on the lattice depth and the pulse duration. The 
lattice depth determines the band structure and is re- 
flected in the terms Cn^q- The pulse duration influences 
the phase evolution of each Bloch state as ^-'^^r^,qTl/h ^ 

The multi-pulse process, which consists of a number of 
single pulses and intervals can be solved as follows. The 
wave function of the condensate after the first pulse ri 
can be derived from Eq. (j4|) as 



l*(n,t)) = ^6(mo,m,ri)e-^^'^'^/^|2mnfeL) 

m 



(5) 



After the first interval r/i and the second pulse r2, the 
population of the \pm) state can be achieved as 



p(2) 



^6(mo,mi,ri)e ^^^/^^(mi, m, r2) 



(6) 



As shown in Eq. (|6]),the population is affected by the two 
pulses, the first one corresponding to 6(mo,mi,ri), and 
the second one corresponding to 6(mi, m, r2). During the 
interval r/i, the phase of the \pmi) state evolves along the 
time as e-^^^""'^^/^, where E^"^^^ = {2mihkL)^/2M = 
Am\ER is the kinetic energy, and Er = {hkLY /2M is 
the single photon recoil energy. The interval r/i pro- 
duces a phase shift e~*^^ ^'^/i/^ and contributes to the 
momentum distribution. 



In the same way, the population of the \pm) state after 
(s + 1) pulses can be achieved as: 



s+1 



(7) 



mi ,m2,--- ,ms i—1 



with mo = 0, and mg+i = m. 

From the analysis above, the momentum distribu- 
tion after a sequence of pulses' scattering is influ- 
enced by not only the lattice pulses with the term 
fl^^^ 6(mi_i, m^, Ti), but also the intervals among the 

pulses as reflected in the term ]^^^-^ e"*^^"^*^'^^^/^. Al- 
though the populations of the momentum states do not 
change during the intervals, the phase-evolution rates of 
the momentum states with different kinetic energies are 
not identical. The phase deviations between the states 
oscillate from to 27r with the interval, and the hetero- 
geneously accumulated phases change the distribution of 
the condensate in the quasi- momentum space. 



III. EXPERIMENTS OF STANDING WAVE 
PULSE SEQUENCES 

We performed the experiments of a condensate in a 
magnetic trap (MT) (see Fig. [U^a)) being scattered by 
a sequence of standing wave pulses (see Fig. mb)). As 
shown in Fig.fljc), after pre-coohng, a cigar shaped ^''Rb 
condensate of 2 x 10^ atoms in 5^5'i/2 \F = 2,Mf = 2) 
state was achieved by the radio frequency (RF) cooling in 
the magnetic trap, of which the axial frequency is 20 Hz 
and the radial frequency is 220 Hz [l5|, [l6|. A pair of 
counter-propagating laser beams, of which the durations 
were controlled by an acousto-optical modulator, and the 
amplitudes were adjusted by the injection current of a 
tapered-amplifier, were applied to the condensate along 
the axial direction. The linear polarized incident light at 
the wavelength = 852 nm was focused with a waist of 
110 /im to cover the condensate. The trap depth, which 
was calibrated by Kapitza-Dirac scattering experimen- 
tally, reached 120£^i^, corresponding to the light power 
of 320 mW. The incident light and the magnetic trap 
were simultaneously shut after the BEC-light interaction. 
After 30 ms free falling and ballistic expansion, the mo- 
mentum distribution of the condensate was pictured by 
absorption imaging. 

The lattice in our experiments is quite deep, so we 
concentrate on the short-pulse diffractions to avoid the 
de-coherence and heating effects of long pulses relevant 
for Bragg scattering. However, for more flexible momen- 
tum manipulation, our pulses are not so short as the 
Raman-Nath pulses [13] used in previous works. 

A brief instruction to the Raman-Nath regime is given 
in the following for comparison. In the scattering process, 
the evolution during the free evolution intervals is ana- 
lyzed as in previous section, while the effect of lattice with 
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FIG. 1: (a) A pair of counter-propagating light beams are 
applied to a condensate in magnetic trap, (b) The scattering 
process consists of a sequence of standing wave pulses, which 
contained s intervals with widths Tfi (i=l,2,...s) and s+l light 
pulses with durations n (i=l,2,...s + 1). The incident light's 
wavelength is 852 nm and its maximum intensity can reach 
120Er. (c) The procedure for the experiments is shown. The 
condensate is exposed to a sequence of standing wave pulses 
and then released from the magnetic trap. The absorption 
images of the condensate can be observed after the free falling. 



adequately short duration r can be analytically solved by 
the Schrddinger equation id\'^{t)) /dt = H\^{t))^ after 
omitting the atomic kinetic energy term /2M in the 
Hamiltonian. This approximation can be made while the 
displacement of the scattered atoms during the interac- 
tion time is much smaller than the spatial period of the 
standing wave. Equivalent ly, the standing wave duration 
r and the single photon recoil frequency uOr = hk\/2M 
have to fulfill r <C l/uor- The pulse is able to split a 
stationary condensate into components with symmetri- 
cal momenta Pn = 2hkL{n = 0, ±1, ±2, ...), with corre- 
sponding populations Pn = J^{Uor/2h)^ where Jn{z) are 
Bessel functions of the first kind. 

First we demonstrate a one-pulse scattering experi- 
ment. A condensate is exposed to a standing wave pulse 
with depth 100 Er and duration varying from to 30/is. 
The relative populations of the condensates with the mo- 
menta Ohk^ ±2fi/c, ±4fi/c and ±6fi/c, corresponding to 
Fig.[2fa), (b), (c) and (d) respectively, are measured and 
theoretically analyzed. In addition, the theoretical anal- 
ysis with the Raman-Nath approximation, is also shown 
in the figure for comparison. It can be seen that within 
3/is the theoretical analysis with the Raman-Nath ap- 
proximation (blue solid line) is close to the experimental 
results (black dots), and so is the theoretical analysis with 
band structure theory (red dashed line). When the pulse 
duration exceeds 3/i5, the analysis with the Raman-Nath 
approximation gradually goes far away from the experi- 
mental results, while the numerical simulation with band 
structure theory still agrees with the experimental results 
along the whole time scale. As shown in Fig. [21 the prob- 
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FIG. 2: Single pulse scattering of condensate: The black dots 
represent the experiment results. The blue solid line is the 
theoretical analysis with Raman-Nath approximation. The 
red dashed line is the numerical simulation with band struc- 
ture theory. Figure (a), (b),(c) and (d) correspond to the rel- 
ative populations of the condensates with the momenta Ohk, 
±2hk, ±4hk and ±6hk respectively. 



ability of each momentum state oscillates with the pulse 
duration as described by the band structure theory. It is 
clear that, in the single pulse scattering process, the band 
structure theory works well not only for the short pulse 
but also for the longer pulse, because the atomic motion 
has been taken into account. So the atomic diffraction by 
a single standing wave pulse can be predicted in a wider 
range of pulse duration with the band structure theory. 

Then we increase the number of pulses in the experi- 
ments to explore the extra factors influencing the momen- 
tum distributions. Two groups of experiments are carried 
out, where one consists in two two-pulse sequences and 
the other uses a train of three pulses or a train of four 
pulses. In every sequence, all the pulses are the same and 
all the intervals are identical to make the experiments 
more convenient to carry out. For further comparison 
between the band structure theory and the analysis in 
the Raman-Nath regime, every single pulse is made short 
enough for the Raman-Nath approximation. 

Two experiments of two-pulse scattering are demon- 
strated in Fig. [3l in which the relative populations of the 
stationary condensate Ponk versus the varied intervals r 
is shown. The parameters of the scattering pulses used in 
different sequences are chosen to be of the same products 
of the lattice depth and the pulse duration, so that each 
pulse affects the condensate equivalently. As shown in 
the figure, the intervals actually affect the final momen- 
tum distribution, and the theoretical analysis with the 
band structure theory and Raman-Nath approximation 
both picture well the evolution of the atomic distribu- 
tions versus the interval between the two pulses. The 
results of two-pulse scattering can be explained as the 



4 



0.° 



I 



U ^ \ \ \ \ \ \ U 

10 20 30 40 50 60 70 80 

interval time 



ex? 



"L \ \ \ ^ X X ^ J 

10 20 30 40 so 60 70 80 

interval time x ^ 

FIG. 3: Two-pulse scattering of the condensate: The relative 
populations of the stationary condensate P^^nk versus the var- 
ied intervals r is shown. The parameters of the experiments 
are described above each figure. The black dots are the ex- 
perimental results. The blue dashed line is the analysis with 
Raman-Nath approximation. The magenta dotted line is the 
numerical simulation with band structure theory. The red 
solid line is a numerical simulation taking into account the 
momentum dispersion. 



fact that since the phase shift accumulated during the 
interval varies harmonically from to 27r, the probability 
of the stationary condensate oscillates between the mini- 
mum and the maximum. When the phase shift is 27r with 
the interval T:h/2ER (around SO/is), the wave function is 
little affected by the interval and the two pulses diffract 
the condensate as one combined pulse to make the prob- 
ability Pohk the minimum. While the phase shift is tt 
with the interval ttH/AEr (about 40/is), the second pulse 
produces an effect opposite of the first one and diffract 
the non-stationary components of the condensate back 
to the stationary one and make the probability Pohk the 
maximum. 

It also can be acquainted from Fig. [3] that the numeri- 



cal simulation with band structure theory is much closer 
to the experimental results than the analytical solution 
with Raman-Nath approximation. It is conjectured that 
the phase evolution during the scattering process makes 
the difference. The phase evolution in the scattering pro- 
cess is neglected in the Raman-Nath approximation, but 
not in the numerical simulation with band structure the- 
ory. Although the duration of the scattering process is 
short, the phase shifts in the scattering process still in- 
crease. The phase shift in the scattering process needs 
to be taken into account and influences the final momen- 
tum distribution. As a result, the longer the scattering 
pulse is, the larger the difference is. Although the maxi- 
mum of the probability Ponk corresponds to the interval 
nh/AEn^ the two-pulse experiments in Fig. [3] can clearly 
show that the longer pulse leads to the larger difference. 
In Fig. [3] (a), the pulse duration is 6/i5, the probabil- 
ity Pohk reaches the maximum with the interval 34/is. In 
Fig.[3](b), the pulse duration is 2/is, the probability Ponk 
gets to the top with the interval 38/is. 

Nevertheless, there is still some obvious deviation be- 
tween the simulation and the experimental results. It is 
observed that the momentum width has been expanded 
after the former pulse, because of the s-wave scattering 
between the different momentum states. Consequently, 
this dispersion process is approximated to an initial mo- 
mentum width of ~ O.lhkL on average to optimize the 
numerical simulation. Unlike the analysis without mo- 
mentum width, phase evolution is different for different 
initial momenta and results in a phase dispersion. The 
quasi modes obtained at the end of the diffraction pro- 
cess result from the linear superposition of final states 
obtained after time evolutions of the different momenta 
populated the initial BEC. It can be seen from Fig. [3] 
that the approximation is effective. 

As discussed in [18], the maximum of the probabil- 
ity Pohk will never reach 1 thanks to the imperfect opti- 
cal lattice. In our case, the momentum expansion is an 
explanation of the similar situation as shown in Fig. [3l 
Since the momentum width is considered, the phase shift 
is populated around tt with a width, instead of a definite 
TT, with the interval ttH/^Er. In other words, there is no 
interval that accumulates a phase shift exactly equal to 
TT, so with any interval, the second standing wave pulse 
is not able to diffract all the condensates back to the 
stationary part. 

The experiments of one three-pulse scattering and one 
four-pulse scattering are demonstrated in Fig. 21 where 
the relative populations of the stationary condensate 
Pohk versus the varied total intervals ^ Tfi is shown. In 
each experiment, the total interaction intensity ^ Uqt 
is the same, with different number of the pulses. These 
two experiments show that even though the total interac- 
tions and intervals are the same, the different processes of 
phase accumulations in the two kinds of pulse sequences 
result in distinct momentum distributions. We directly 
apply the band structure theory with momentum dis- 
persion to analyze the experiments in Fig. 31 and the 
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FIG. 4: Three-pulse and four-pulse scattering of the conden- 
sate: The relative populations of the stationary condensate 
Pohk versus the varied total intervals X^r/^ is shown. The 
parameters of the experiments are described above each fig- 
ure. The black dots are the experimental results. The red 
solid line is a numerical simulation taking into account the 
momentum dispersion. 



corrected simulations agree with the experiments quite 
well. 



IV. MANIPULATE THE MOMENTUM STATES 
AS DESIGN 

The experiments and the numerical simulations above 
have shown the possibility and feasibility of the manip- 
ulation of a condensate's momentum states. We man- 
age to design several two-pulse sequences to achieve 
high contrast momentum states such as |±2fi/c), |±4fi/c) 
and \^^hk), which may be useful in atomic interferome- 
try [l9|,[20|. For each state, we apply two totally different 
two-pulse sequences to show the flexibility of the method. 
The general method to achieve the target states is to find 
out the condition of the minimum of the square deviation 

= E (^m - ^m) , where is the probability of 

m— — oo 

\2mhkL) in the goal state, and Pm is that generated by 
the sequence. A second method, as the target is to obtain 
the highest population of some certain momentum state, 
consists in scanning the whole set of initial conditions 
and choose the one corresponding to the maximum value 
of the desired population. We apply the two methods 
above separately and obtain the same pulse sequences. 
As shown in Fig. [5l the experimental results (the black 
round dots) agree well with the expectations of the de- 
signs (the blue diamond dots), whether the pulses are in 
the Raman-Nath regime (see Fig. Efb)) or not (see oth- 
ers in Fig. [5]). When the momentum dispersion is being 
considered, the expected momentum distributions (the 
red square dots) get closer to the experiments, where the 
figures only display the relative populations of the tar- 
get states and omit the others for the figures being more 
clear. 

It can be seen from Fig. [5] that the momentum width 
correction can improve the precision of the prediction 



momentum p(2%k) 

FIG. 5: Experimental realization of designed momentum 
states. The expected momentum state is ±.2hkL{{^) and (b)), 
±4/iA;l((c) and (d)) and ±6/i/cL((e) and (f)). The pulse se- 
quences are shown above each figure. The black round dots 
are experimental results. The blue diamond dots are the ex- 
pectations based on the design. The red square dots are the 
modified design with momentum width, which agree with the 
experiments better. 



with our method. The average relative deviation between 
the experimental results and the expected values without 
including the momentum width is 25.03%, while the de- 
viation is decreased to 13.15% with the correction. 

An asymmetry of the momenta can be observed in 
Fig. [5l and it may ascribe to the following factors. Be- 
sides the measurement error, there is an imperfection of 
the standing wave, brought forth by the unbalanced in- 
tensity of the laser beams. External field (such as the 
magnetic trap) fluctuations during the scattering process 
may also affect the momentum distribution. 



V. DISCUSSION AND CONCLUSION 

The band structure theory is a global method to deal 
with the standing wave scattering a condensate, while 
Bragg and Raman-Nath scattering are two special situa- 
tions which can be analytically solved with their respec- 
tive approximations. In the Bragg regime, the poten- 
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tial height introduced by the standing wave is restrained 
below 4:Eji and that leads to the difficulty of generat- 
ing higher order momentum states. In the Raman-Nath 
regime, the intensity of the standing wave is not limited 
so that higher order momentum states can be generated 
symmetrically [2l|, [22| . However the pulse duration has 
to be short enough to neglect the atomic motion, so the 
momentum states can not be predicted in this regime if 
the pulse duration is slightly longer. In our work, the 
scattering can be well explained and numerically ana- 
lyzed within a much wider range of pulse intensity and 
duration. So it is natural that more momentum states 
can be generated. 

In our article, we compared the scattering by one pulse 
and that by a train of pulses. Some valuable states, such 
as |±2fi/c), |±4fi/c) and \^6hk) states with high contrast, 
can not be realized by the single pulse scattering, while 
they can be realized by a sequence of standing wave 
pulses. A sequence of lattice pulses is a more effective 
and flexible tool for momentum manipulation. It can 
generate much more useful momentum states than the 
ones demonstrated in our work. In the future, more pa- 
rameters could be changed for the design to obtain better 
results. 

Even though the numerical simulation is corrected to 
take into account the momentum dispersion, some devi- 
ations between the experiments and the simulation still 
exist. The inaccuracy of the lattice-depth calibration, 
which is 5% at least, is one of the reasons. The phase shift 



introduced by the magnetic trap is another one, while the 
influence is estimated to be within 0.03%, which is below 
the experimental uncertainty. The heating and momenta 
exchange during the s-wave scattering may also lead to 
some differences, which need further study. 

In conclusion, we developed a method for more flexi- 
ble manipulation of the condensate's momentum states, 
where the momentum states can be controlled by stand- 
ing wave pulses in a wider range of pulse intensity, du- 
ration or quantity. The experiments of a condensate be- 
ing scattered by a sequence of standing wave pulses are 
demonstrated. A global theory, treating the standing 
wave as an optical lattice, is applied to explain the ex- 
periments. With this theory, we are able to design pulses 
sequences for realizing states such as \±2hk), \±4:hk) and 
|±6fi/c), and experimentally realize them, which may be 
applied in atomic interferometry to improve the measure- 
ment precision. 
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